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1 Introduction 
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The full next-to-leading order (NLO) corrections to the hadroproduction of 
heavy flavors have been completed in 1988 [1, 2]. They have raised the leading 
■ order (LO) estimates [3] but were still below the experimental results (see e.g. 

\ [4]). In a recent analysis theory moved closer to experiment [4]. A large un- 

^sO ■ certainty in the NLO calculation results from the freedom in the choice of the 

CD renormalization and factorization scales. The dependence on the factorization 

and renormalization scales is expected to be greatly reduced at next-to-next- 
to-leading order (NNLO). This reduces the theoretical uncertainty. Further- 
f^S ■ more, one may hope that there is yet better agreement between theory and 

experiment at NNLO. 

In Fig. 1 we show one generic diagram each for the four classes of contribu- 
tions that need to be calculated for the NNLO corrections to hadroproduction 
of heavy flavors. They involve the two-loop contribution (Fig. la), the loop-by- 
pH | loop contribution (Fig. lb), the one- loop gluon emission contribution (Fig. lc) 

and, finally, the two gluon emission contribution (Fig. Id). An interesting sub- 
class of the diagrams in Fig. lc are those diagrams where the outgoing gluon 
is attached directly to the loop. One then has a five-point function which has 



to be calculated up to 0(e 2 ). 

In our work we have concentrated on the loop-by-loop contributions exem- 
plified by Fig. lb. Specifically, working in the framework of dimensional reg- 
ularization, we have calculated 0(e 2 ) results for all scalar one-loop one-, two- 
and three-point integrals that are needed in the calculation of hadronic heavy 
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Fig. 1. Exemplary diagrams for the NNLO calculation of heavy hadron production 

flavour production. Work on the relevant four-point integrals is in progress. 
The integrations were generally done by writing down the Feynman parame- 
ter representation for the corresponding integrals, integrating over Feynman 
parameters up to the last remaining integral, expanding the integrand of the 
last remaining parametric integral in terms of e and doing the last parametric 
integration on the coefficients of the expansion. Because the one-loop inte- 
grals exhibit infrared (IR)/collinear (M) singularities up to 0(e~ 2 ) one needs 
to know the one- loop integrals up to 0(e 2 ) because the one- loop contribu- 
tions appear in product form in the loop-by-loop contributions. It is clear 
that the spin algebra and the calculation of tensor integrals in the one-loop 
contributions also have to be done up to 0{e 2 ). This work is in progress. 

Due to lack of space we can only present a few exemplary results. Since 
the four-point functions are the most difficult we concentrate on them. 



2 Four-point functions 

As a sample calculation we discuss the (0, 0, m, 0)-box with one massive prop- 
agator depicted in Fig. 2. As explained before one needs to calculate each 




Fig. 2. One-loop box with one internal massive propagator 



one-loop integral up to 0(e 2 ) in order to obtain the finite terms in the loop- 
by-loop contributions. The box integral Fig. 2 is represented by the integral 
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D(-p2,P4,P3,0,0,m,0) = 



* J (27T)" 



(q 2 ){q -P2) 2 [{q -Pi +P4) 2 - m 2 ](q - p 2 + p 4 + P3) 1 



where pi,P2, Pz and p 4 are external momenta with p\ = p\ = 0, p\ = p\ = m 2 
and n = 4 — 2e is the dimension of space-time. 

The e~ 2 , e _1 and e° coefficients have been known for some time [1, 2] 
and will not be listed here. We define Mandelstam-type variables by s = 
(pi+p 2 ) 2 ,t=( Pl -p 3 ) 2 -m 2 and/3= {I - Am 2 / s) 1 / 2 , x = (1 - +/?). 

For the real part of the 0(e) term we obtain: 
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6 In 3 A + 20 In 3 + In 3 x + 6 In 2 x In ' (s + 2t ~ s;3) 



_ 12 In* In 2 zls^lM ) + 8 ^ -(*+2t-sp) 

2 m 2 I 2 m 2 

-12 1n 2 ^| f2 ln(-l - A) + 3 Ins - In zLs^±3 + 4 i n *+2*+»fl 

V v vn z > 2m z 2m A 

_ 3 l n 2 _2, f 6 In -4 + 3 In a: + 2 In -('+";" + 4 In a+ "+^ > > + f 2 In + 



5 : h" 4: 111 „ 5- 

24 ln ^| + 24 ln ^ ln ^ + 48 Li 3 {^) + 48 Li 3 (l + + 24 Li 3 (-x) - 
24 Li 3 C 2m ' t $v 1+(i) ) 48 L, 3 (^±^) + 24 L* 3 (^+*^) 

- 24Li 8 (^^) ~ tfLiafjis^) + 24Li 3 (-^^i±g) 
-24^ 3 (^^)-24^ 3 ( 2 ^±^)-6 ln^| (4+2 ln 2 (-l-^)-ln 2 x 
4 In a; ln -(«+ 2 «-^> + 6 In 2 z(£+g^) 

2 2 m J 

- 4 ln(-l - (in* + In + In + 8 Li 2 (l + 
+ 8L. 2 (i™) + 4L i2 (^i±g J ) -4^ 2 ( ^)^ ) 
+8^(OTW)- 12 C(2))-12 lna:C(2)+48 ln -(^-^) c(2) _ 6 m ^ ( 2 + 
41n 2 ^-ln 2 x-2 1nxln^|^-41n^(-f + lnx + ln^±i^ + 
2 ^ - 4 «a(±^) - 4 Li 2 (^if^) + 2 C (2)) - 36 C (3) , 
where C e (m 2 ) ^ ^ 

One also needs the imaginary part of the (0, 0, m, 0)-box since the to- 
tal contributions from the loop-by-loop contribution contains also imaginary 
parts via \A\ 2 = (ReA) 2 + (ImA) 2 . Note , however, that the imaginary part is 
only needed up to 0(e) since the IR/M singularities in the imaginary parts of 
the one-loop contributions are of 0(^) only. For the 0(e) absorptive (imagi- 
nary) part we obtain: 



3 In 2 A -4 ln A m^|-8 In 2 =\-2 ln lnx + 4 ln =\ lnx 
3 ln 2 x - 4 ln A In ^S^l3 _ 4 lnx ln zis+^lR + 4 In 2 ziS^lB 

ra z 2m z 2m z 2m z 
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4 1n^ln^^ + 8 1n^ln^^+4 1 na; ln^^ + 8L i2 (^)) + 
8£i2( sJfe? ) + 8C(2) • 



The e 2 -results for the (0, 0, m, 0)-box are too lengthy to fit into this report. 
They will be presented in a forthcoming publication [5]. A new feature of 
the e 2 - contributions is that the result can no longer be expressed in terms 
of logarithms and polylogarithms . They involve more general functions - 
the multiple polylogarithms introduced by Goncharov in 1998 [6]. A multiple 
polylogarithm is represented by 

x 1 x 2 ...x k _ TOl -l , 

at \ dt 



Li mk ,...,mi(Xk, -,xi) = / — 



o 



t J X2X 3 ...X k - t 



dt X" 12 " 1 dt (dt \ m *- x dt 



o o ... o — o 



t ) X 3 ...Xk—t \t J l—t/ 

where the iterated integrals are defined by 

a a t n t 2 



/dt dt f dt n f dt n f 
o ... o = / / x ... x / 
a n — t a\—t J a n — t n J a n _i — t n -\ J 



dt 



ai - ti 
o 



Besides the scalar (0, 0, m, 0)-box one also needs to calculate the scalar 
(0,0, m,m) and (m, m, 0, m)-boxes. Work is in progress on the calculation 
of these boxes. 



3 Summary 

We have reported on the results of an ongoing calculation of the e-expansion 
of the scalar one- loop integrals up to 0(e 2 ) that are needed for the NNLO 
calculation of hadronic heavy hadron production. In order to arrive at the 
full amplitude structure of the one-loop contributions one still has to include 
the £-dependence resulting from the Passarino-Veltman decomposition of the 
tensor integrals, and the ^-dependence of the spin algebra calculation. Putting 
all these pieces together one might optimistically say, when considering the 
four classes of diagrams Fig. 1, that the present calculation constitutes one- 
fourth of the full NNLO calculation of hadronic heavy hadron production. 
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